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An  Incremental  Form  of  the  Single-Integral 
Nonlinear  Viscoelastic  Theory  for  Elastic- 
Plastic-Creep  Finite  Element  Analysis 

D.  a SANDERS  W.  E.  HAISLER 


ABSTRACT 

single-integral  nonlinear  viscoelastic  theory 
Is  cast  Into  a generalized  Incremental  form  for  use  In 
predicting  creep  response  of  materials.  The  Incre- 
mental creep  constitutive  model  Is  tailored  for  use  In 
an  Incremental  finite  element  program  wherein  the 
equations  of  motion  are  Integrated  numerically  step- 
by-step.  Recursive  relations  for  material  memory  pa- 
rameters are  developed  which  allow  the  retention  of 
all  past  history  prior  to  the  current  time  step  and 
require  Integration  over  the  current  time  step  only. 

The  creep  model  parameters  required  In  the  theory  are 
obtainable  from  a standard  creep  and  recovery  test. 

The  Incorporation  of  the  model  within  an  Incremental 
elastic-plastic  finite  element  formulation  Is  outlined. 

INTRODUCTION 

In  a majority  of  the  finite  element  programs  In 
use  today  for  analyzing  metal  components,  the  tradi- 
tional approach  of  separating  the  total  strain  Into 
elastic,  plastic,  and  creep  components  Is  used  (1,2,3). 
The  most  widely  accepted  means  of  predicting  the  creep 
strain  In  this  context  is  the  phenomenological  creep 
theory  (4,5).  As  an  alternative  to  this  approach,  a 
single  Integral  nonlinear  viscoelastic  theory  has 
been  obtained  by  Schapery  (6)  from  thermodynamic  con- 
siderations. It  has  a form  similar  to  the  Boltzman 
Integral  used  In  linear  viscoelasticity,  but  the 
physical  time  Is  replaced  by  a so  called  reduced  time; 
and  several  other  material  functions  enter  Into  the 
theory.  In  theory,  this  representation  of  the  Inelastic 
strain  components  should  provide,  within  the  same 
constitutive  theory,  the  capability  of  analyzing  both 
polymeric  and  metallic  materials. 

Two  constitutive  theories  similar  In  appearance 
to  the  single  Integral  theory  of  Schapery's  are  given 


by  Rashid  (7)  and  Valanls  (8).  Rashid  uses  a modified 
superposition  principle  with  the  compliance,  which  Is 
a function  of  stress  and  time,  transformed  such  that 
It  Is  a function  of  a single  parameter  called  reduced 
time.  Valanls's  theory  (8)  Is  derived  from  thermody- 
namic considerations  with  the  resulting  single  Integral 
equation  appearing  In  terms  of  the  strain  history  In- 
stead of  the  stress  history  as  In  reference  (6).  This 
theory  Is  very  similar  to  Schapery's  theory  (9)  wherein 
strain  Is  employed  as  the  Independent  variable  Instead 
of  stress.  The  functional  difference’  In  references  (8) 
and  (9)  Is  due  to  their  respective  definitions  of 
reduced  time.  The  other  material  functions  entering 
Into  references  (6,9)  are  not  present  In  Rashid's  or 
Valanls's  theory. 

Application  of  the  single  Integral  theory  has 
been  made  by  Valanls  (10,11)  for  the  viscoplastic 
response  of  several  metals.  It  has  also  been  shown  to 
give  accurate  results  for  several  polymeric  composites 
In  creep,  creep-recovery  and  multlstep  Input  tests  on 
uniaxial  specimens  (12,13). 

A potential  advantage  of  the  single  Integral 
approach  Is  that  It  easily  handles  multlstep  Input 
histories  Including  stress  reversal  and  creep-recovery. 
This  Is  not  true  of  the  phenomenological  theory  of 
creep  for  metals  based  on  strain  hardening  which  must 
resort  to  special  procedures  to  handle  the  stress 
reversal  situation  as  outlined  In  (14).  These  proce- 
dures have  In  some  cases  resulted  In  rather  poor  results 
as  Indicated  In  reference  (14).  Further,  the  phenome- 
nological theory  has  no  predictive  capability  for  the 
recovery  response  to  stress  removal . Thus,  the  single  Inte- 
gral constitutive  theory  may  offer  additional  predic- 
tive capabilities  over  the  phenomenological  theory  of 
creep. 

This  paper  casts  the  single  Integral  theory  Into 
a form  suitable  for  use  within  an  Incremental  finite 
element  program.  Since  the  equations  of  motion  are 
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Integrated  numerically  step-by-step,  the  time  deriv- 
ative of  the  single  integral  equation  must  be  ob- 
tained. In  this  form,  the  creep  strain  Increment  can 
be  determined  for  the  current  time  interval  based  on 
the  creep  strain  rate  at  the  beginning  cf  the  step. 
Computationally,  the  relationships  derived  for  the 
rate  form  of  the  single  integral  equation  are  partic- 
ularly appealing  since  recursive  relationships  are 
developed  which  allow  for  the  retention  of  ail  past 
history  prior  to  the  current  time  step  and  require 
integration  over  the  current  time  step  only. 

The  incorporation  of  the  model  within  an  incre- 
mental elastic-plastic  finite  element  formulation  is 
outlined.  The  time- Independent  Incremental  plasticity 
relations  with  a von  Mises  yield  criterion  and  a 
kinematic  hardening  rule  are  employed  to  describe  the 
plastic  component  of  strain. 

A graphical  procedure  for  determining  the  mate- 
rial functions  from  either  creep  data  or  creep  equa- 
tions Is  presented.  It  is  shown  that  this  procedure 
is  only  applicable  to  materials  whose  elastic  strain 
change  on  unloading  Is  equal  to  that  on  loading. 

Numerical  comparison  of  the  single  integral 
theory  to  that  of  phenomenological  theory  of  creep 
are  given  for  creep,  creep  and  recovery,  multistep 
and  reverse  loading  histories  for  a type  304  stain- 
less steel.  The  numerical  results  were  obtained  using 
a computer  program  NONVIS  developed  by  the  authors 
which  Implements  the  Integration  procedure  developed 
below  for  the  single  Integral  creep  equation. 

It  is  suggested  that  the  rate  form  of  the  single- 
integral equations  appearing  In  (7,8,9)  could  be  Inte- 
grated by  employing  a procedure  similar  to  that  given 
In  this  paper  for  the  single  Integral  constitutive 
theory  of  reference  (6),  provided  the  kernel  function 
Is  a combination  of  constants,  linear,  and  exponential 
terms. 

RATE  FORM  OF  THE  SINGLE  INTEGRAL  EQUATION 

The  single  integral  nonlinear  viscoelastic  con- 
stitutive equation  derived  In  reference  (6)  for  an 
Initially  Isotropic  material  Is  written  for  Isothermal 
conditions.  Incompressible  creep  strains,  and  Infin- 
itesimal strains  as. 
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and  A'1',  Bv  ' and  rP  are  constants.  The  functions 
ag.aQ  and  3C(j  are  material  properties  and  are  In 
general  functions  of  the  three  stress  Invariants. 

The  quantity  * Is  a reduced  time  variable  while  ca(j 
are  engineering  strain  components  and  oC(j  are  compo- 
nents of  the  Cauchy  stress  tensor.  The  function  aJ(p) 
Is  the  transient  component  of  the  so  called  linear 
viscoelastic  shear  creep  compliance.  The  first  term 
In  equation  (1)  represents  the  instantaneous  elastic 
response  of  the  material. 


To  cast  equation  (1)  Into  a form  suitable  for  use 
In  an  incremental  finite  element  analysis.  Its  time 
derivative  Is  taken,  and  equation  (2a)  Is  used  to 
obtain, 
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Equation  (3)  Is  rewritten  as, 
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A time  marching  procedure  for  evaluating  Xcd  and  Ycd 
Is  formulated  below. 


Evaluation  of  XC()  and  Y ^ 

The  evaluation  of  the  two  Integrals  Xcd  and  Ycd 
given  by  equations  (5)  and  (6)  represents  a key  element 
In  evaluating  the  creep  strain  rate.  Once  the  value  of 
these  Integrals  are  known,  the  creep  strain  rate  Is 
evaluated  from  the  sum  of  the  second  and  third  terms 
of  equation  (4)  as: 
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Two  Important  features  In  the  evaluation  of  Xcd  and  Ycd 
are  the  determination  of  recursive  formulas  which 
make  use  of  values  from  the  previous  time  step  and 
require  the  evaluation  of  the  two  Integrals  over  the 
current  time  step  only.  The  formulation  presented 
here  Is  similar  to  the  Integration  technique  employed 
by  Zak  (15)  and  Bills,  et  al.  (16)  for  a thermorheo- 
loglcal  simple  material  due  to  a variable  temperature 
field. 

Assume  the  time  from  0 to  t Is  subdivided  Into 
N subintervals  such  that  to*0,t]«to*Ari , ...» 
tN-tN.i+atN-t.  For  ocd(x.tN)  and  Erdk.tN).  write 
ocdN  and  GCdu  and  assume  that  ocdwt^^N)  * 0 *nd 
GcdjJfO.tfl)  = 0.  Then  Xcd  and  Ycd  can  be  written  as 
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Evaluation  of  X/-^.  Using  equation  (2c)  In 
equation  (9),  one  obtains 
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By  evaluating  the  first  Integral  In  equation  (11), 
expanding  the  second  and  third  as  the  sum  of  two 
Integrals  and  approximating  the  time  rate  of  change 
of  G^,  within  the  time  Interval  t^.^^t^  as 
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then  equation  (11)  can  be  written  as 
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Define  the  following  quantities. 
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With  the  use  of  equation  (16),  equation  (13)  Is 
rewritten  as 
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Recursive  relationships  for  equations  (23)  and  (24) 
are  developed  In  the  Appendix  and  are  given  below  as: 
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where  Cl^Jj^,  C2^dN-l*  J1$T|  *nd  known  from 


the  previous  time  step.  Thus,  equation, (19)  can  be 
evaluated  by  Integrating  Jljjj)  and  J2^'def1ned  by 
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equations  (17)  and  (18),  respectively,  over  the  current 
time  step  only.  All  history  effects  are  contained  In 
lCdN  defined  by  equation  (20).  In  the  Appendix, 
equations  (17)  and  (18)  are  numerically  Integrated  using 
Gaussian  quadrature.  Equation  (17)  becomes: 
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Evaluation  of  tho  reduced  time  Integral  given  In 
equation  (2a)  Is  also  presented  In  the  Appendix. 


Evaluation  of  The  Integration  of  YC(jn  Is 

carried  out  In  a similar  manner  to  that  used  for  Xcd 

and  Is  outlined  briefly  here.  rf(L  was  given  by 
equation  (10)  as,  N 
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Expanding  oquat  on  (10)  as  the  siaa  of  two  Integrals, 
using  equation  (31)  and  approximating  the  time  rate  of 
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Mote  that  Jl^  Is  the  same  quantity  that  appears  In 
the  evaluation  of  Xcd.  Using  equation  (33)  and  (34), 

equation  (32)  Is  rewritten  as. 
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A recursive  relationship  for  equation  (37)  Is  developed 
In  the  Appendix  and  Is  given  below  as: 
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*•**»“•  F^  ^ and  Jlj^j  are  known  from  the  previous 

time  step.  Thus  Ycd^,  In  equation  (37),  can  be 
evaluated  by  Integrating  Jl£r*  given  In  equation  (17) 


over  the  current  time  step  only  using  numerical 
quadrature  as  outlined  In  the  previous  section.  All 
history  effects  are  represented  by  given  by 
equation  (34).  " 


INCREMENTAL  EQUATION  OF  EQUILIBRIUM 


In  this  section,  an  Incremental  finite  element 
formulation  for  the  Isothermal,  elastic-plastic-creep- 
large  strain  problem  Is  developed.  The  formulation 
1$  suitable  for  use  with  any  creep  constitutive  theory 
which  will  predict  the  Incremental  creep  strain  com- 
ponent based  on  the  creep  strain  rate  at  the  beginning 
of  the  step.  Although  the  creep  and  plasticity  rela- 
tions are  for  small  strains,  the  formulation  Is  devel- 
oped for  large  strain  kinematics.  One  may  begin  with 
the  equations  of  equilibrium  written  In  terms  of  the 
second  Plola-KIrchhoff  stress: 
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where  aj  and  u,  are  Langranglan  coordinates  and  dis- 
placements, respectively,  Sjk  Is  the  2nd  Plola- 
Xlrchhoff  stress  tensor,  p0  Is  undeformed  density  and 
F0j  Is  the  body  force  per  unit  undefonned  volume  per 
unit  mass.  Applying  the  virtual  work  principle  at 
time  t*at  yields 
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•ltd  SEfi  Ms  the  variation  of  the  Green- La grange 
strains  at  time  t*At  and  T?+itare  surface  tractions 
at  time  t+at  applied  to  thl  deformed  surface  S. 
Equation  (40)  may  be  put  Into  Incremental  form  by 
writing 
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where  sfj  and  eJj  are  stresses  and  strains  at  time  t 
and  aSjj  and  aE^  are  Increments  of  stress  and  strain, 

respectively.  The  strain  Increment  may  be  decomposed 
Into  components  which  are  linear  and  nonlinear  In  the 
displacement  Increments 
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Substituting  equations  (42)  and  (43)  Into  (40)  yields 
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The  stress  Increment  may  be  decomposed  Into  two  com- 
ponents, one  which  Is  dependent  upon  total  strain  and 
one  which  Is  Independent  of  strain  (l.e.,  creep, 
thermal,  etc.): 
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DlJkl  1*  the  usual  effective  tangent  modulus  ai 
aPjj  Is  i stress  Increment  due  to  strain  Independent 
phenomena  (as  Is  usually  assumed  In  creep).  Substi- 
tuting equation  (46)  Into  equation  (45),  making  use 
equation  (43),  and  neglecting  terms  which  would  be 
nonlinear  In  displacement  Increments,  yields  the 
following: 
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The  term  aPji  may  be  Interpreted  as  the  change  in 
stress  required  to  account  for  the  creep  and  thermal 
strains.  Equation  (47)  takes  on  the  following  form 
when  put  Into  matrix  form 

[M]  (q4*44)  ♦ {[K4]  ♦ [KjL]){aq)  - (R4*44)  - <F4) 

(48) 

where  [M]  Is  the  mass  matrix,  [K4]  and  [K4  ] are 
"linear"  and  nonlinear"  stlffnest  matrices,  (Rt*At>  i$ 
a vector  of  forces  due  to  externally  applied  loads, 
IF*)  Is  a vector  of  forces  due  to  Internal  stress,  and 
(Aq)  Is  the  Increment  of  the  nodal  displacements. 
Complete  details  of  the  derivation  of  the  quantities 
In  equation  (48),  without  creep,  may  be  found  In 
reference  (17). 

The  determination  of  Dniti  and  aPji  depends  on 
the  specific  hardening  rule  chosen.  For  simplicity 
of  presentation,  the  Isothermal  model  Is  chosen  and 
It  Is  assumed  that  the  yield  function  can  be  expressed 
by 
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For  kinematic  hardening, 
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where  S<j  and  are  devlatorlc  components  of  t 
Plola-Klfchhoff  stress  and  yield  surface  center, 
respectively  and  K Is  a constant. 

The  associated  flow  rule  Is  given  by 
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where  X Is  a scalar  to  be  determined.  Using  the 
Ziegler  modification  of  the  Prager  kinematic  work- 
hardening rule,  it  Is  assumed  that  the  projection  of 
the  Increment  of  stress  onto  the  normal  to  the  yield 
surface  Is  a scalar  multiple  of  the  dot  product  of 
the  plastic  strain  and  yield  surface  .normal : 
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where  C Is  a scalar  (hardening  modulus)  to  be  deter- 
mined from  a uniaxial  stress-strain  curve.  It  can  be 
shown  that  equation  (52)  Is  equivalent  to  the  consist- 
ency condition  for  kinematic  hardening  and  Isothermal 
conditions.  For  small  strains,  a decomposition  of 
strain  Is  assumed: 
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where  E1Jmn  Is  the  elastic  constitutive  tensor  and 
dE^,  are  trie  Incremental  creep  strain  components.  The 
Incremental  creep  strain  components  are  determined 
from  equation  (8)  for  a time  Increment  over  which  It 
Is  assumed  the  creep  strain  components  are  constant. 
Substituting  equation  (51)  Into  equation  (52)  and  (53) 
and  solving  the  resulting  equation  for  x gives 
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Substituting  equation  (54)  Into  equation  (51)  and  the 
results  Into  equation  (53)  and  comparing  the  results 
with  equation  (46).  one  can  show  that  the  Instanta- 
neous modulus  tensor  Is  given  by 
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and  the  change  of  stress  due  to  creep  strains  Is 
given  by 
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Equation  (53)  can  now  be  rewritten  as. 


dS 


1J 


D1jkl<dE,«f<n>  • 


DETERMINATION  OF  MATERIAL  FUNCTIONS  FROM  EXPERIMENTAL 
DATA 

In  this  section,  equation  (1)  Is  specialized  to 
the  case  of  uniaxial  loading;  and  to  materials  whose 
elastic  response  on  loading  and  unloading  In  a creep 
and  recovery  test  are- the  same.  The  required  material 
functions  can  then  be  evaluated  by  a graphical  pro- 
cedure from  the  results  of  creep  tests.  For  a gen- 
eral description  of  the  methods  for  obtaining  the 
material  functions  from  experimental  data,  see  Lou 
and  Schapery  (13). 

Consider  a uniaxial  bar  which  Is  In  a uniform 
uniaxial  state  of  stress  at  a uniform  temperature. 
Considering  only  one  stress  o,  and  one  strain  c.  we 
let  cacu  and  When  equation  (1)  Is  special- 

ized to  this  uniaxial  loading,  the  creep  term  Is 
given  by 


r 


d(g2o) 

AD(e-e)  —r~—  dv 


(58) 


and  the  three  Independent  material  functions  are  given 
by 

do. 


9l(o)  " ' 


(59) 


g2o  ■ 2- 
c aG 


(55) 


and 


ao(o) 


. TO 

’ a^To) 


(64) 


(65) 


(56) 


(57) 


When  equation  (63)  Is  subjected  to  a creep  test, 
the  resulting  equation  can  be  used  to  evaluate  the 
material  functions  g2  and  a„.  Consider  equation  (63) 
subject  to  the  stress  history: 

o(t)  - H(t)  o . (66) 

Substituting  equation  (66)  Into  equation  (63)  gives 

tC  - 92  a0  (J-)o  ; (67) 

O 

and  the  creep  coippl lance  Is  given  by 

4D(t)  - 4 ’ 92  A0(a_)  • (68) 

0 

Taking  the  log  of  both  sides  of  equation  (68)  gives, 

log  aD(t)  ■ log  g2  ♦ log  aO(j-)  . (69) 

0 

The  simple  form  of  equation  (69)  Is  such  that  a simple 
graphical  procedure  can  be  employed  to  evaluate  the 
material  functions.  This  Is  easily  demonstrated  by 
plotting  aD  and  t on  double  - logarithmic  paper. 

Curves  at  different  stress  levels  can  be  superimposed 
by  translating  them  along  the  aD  and  t axes.  If  the 
curves  are  shlfceu  *o  a reference  stress,  or,  the 
amount  of  horizontal  shift  (t)  and  vertical  |h|ftj(AD) 

equals  log  g2  • log  and  log  *g  - log 

respectively.  These  shifts  are  graphically  demon- 
strated In  Figure  (1). 


g2(o)  - Gn  (o)  • Jll  , 
b 


TO 


a0(°) 

‘TO 


and 


(60) 

(61) 


AD(e)  - E A(r)  (l-e*<’/Tr)  ♦ . (62) 

It  Is  shown  In  reference  (12)  that  gj(o)  ■ 1 implies 
that  the  amount  of  elastic  straining  on  unloading  is 
the  same  as  on  loading  and  that  “-•'ab*  Using 
91  (a)  ■ 1 and  <iab"aab<  equation  (58)  is  now  given  by. 


c d(g2u) 

EC  • / 4D(*-e')  -dTZ  dT 


(63)  Fig.  1 Horizontal  and  Vertical  Shift  of  Creep  Compliance 
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The  direction  of  the  shifts  are  as  follows: 


92(«) 

iprsp  * 1 * sh1ft  15  up* 

fcto) 

g ^ | > 1 ■ shift  Is  down, 

and 

• (or) 

a < 1 , shift  Is  to  left. 


• 92(or) 

(73c) 

and 

T1  " Vao(or) 

(73d) 

Rearranging  equation  (73a)  gives 

A/o  - g2 

(74) 

and  from  equation  (2a)  and  (73b) 

r • 1/a 

0 

(75) 

• («r) 

■ » 1 , shift  Is  to  right. 

SOME  NUMERICAL  RESULTS 

This  section  presents  some  maaerlcal  results 
obtained  for  a type  304  stainless  steel  with  the 
computer  program  NONVIS  developed  by  the  authors 
which  Implements  the  procedures  developed  above  for 
Integrating  equation  (1)  for  the  following  stress 
histories:  creep  (figure  2) .creep-recovery  (figure  3), 
multi-step  (figure  4),  and  stress  reversal  (figure  5). 
The  results  obtained  are  compared  with  those  of  the 
phenomenological  creep  theory  based  on  strain  harden- 
ing using  total  creep  strain  (IS). 

An  effective  creep  equation  of  the  form: 

:c  - A(1  -e_rt)  ♦ kt  (70) 

was  used  to  determine  the  material  functions  In  lieu 
of  the  actual  experimental  data.  Two  different  heats 
of  type  304  stainless  steel  are  considered.  The  creep 
strain  parameters  for  heat  no.  9T2796  at  llOOf  are 
taken  as: 


A(S)  • 5.436  x 10*5  51*043 

(71a) 

r(o)  ■ 5.929  x 10'5  exp(0.2029  5) 

(71b) 

k(5)  - 6.73  x 10-9  [s1nh(0.1479  5)]3,0 

(71c) 

with  S(effect1ve  stress)  In  ksl,  and  t In  hour-.  The 
creep  strain  parameters  for  heat  do.  9043813  at  1200f 
are  taken  as 

A(o)  • 2.33  x 10‘6  S3-083 

(72a) 

r(S)  • 1.354  x 10"3  exp(0.129  5) 

(72b) 

and 

k(o)  - 7.91  x 1.0*n  [s1nh(0.1932  S)]4’0  . 

(72c) 

The  vertical  shift  Is  defined  by  the  function 

A(o)  and  the  horizontal  shift  by  the  function  r(S). 
This  can  be  seen  be  equating  the  transient  term  given 
In  equation  (58)  (for  r • 1,  r*  - 0)  and  (70)  from 
which  It  may  be  shown  that 

A - g2o 

(73a) 

♦ ■ rt 

(73b) 

Therefore,  the  actual  shifting  of  the  creep 
compliance  1$  not  necessary  to  determine  the  material 

functions  g2(o)  and  jnryx-  rt  is  eas*  t0  validate  the 
1 0 

function  g2  and  — given  above  by  using  the  graphical 
procedure.  a° 

If  the  steady  state  creep  term  In  equation  (70), 
l.e.,  the  term  linear  In  t.  Is  subjected  to  the  creep 
stress  history  given  by  equation  (66),  one  obtains 


(76) 


Using  equations  (74)  and  (75),  equation  (76)  becomes 

cs  - Art  . (77) 

Equation  (70)  above  gives  for  the  steady  state  creep 
term 


(78) 


However,  from  equations  (77),  (78),  and  (71c)  or  (72c), 
It  Is  seen  that  equations  (77)  and  (78)  are  not  equal. 
Thus  It  Is  concluded  that  a single  reduced  time  cannot 
be  used  to  relate  both  the  transient  and  steady  state 
creep  terms  using  the  creep  equations  given  by  either 
equation  (71)  or  (72).  Hence,  In  the  present  appli- 
cation to  the  type  304  stainless  steels,  the  steady 
state  term  will  be  taken  as  the  second  term  of  equation 
(70),  es  ■ kt. 

figure  2 presents  creep  strain  vs.  time  for  a 
uniaxial  stainless  steel  specimen  (heat  No.  9T2796) 
for  a creep  test.  The  results  of  both  the  nonlinear 
viscoelastic  and  the  phenomenological  creep  theory  are 
nearly  Identical  as  Is  ecpected.  Thus,  only  the  non- 
linear viscoelastic  theory  Is  presented. 

figure  3 presents  the  creep  strain  response  of  a 
type  304  stainless  steel  (heat  No.  8043813)  uniaxial 
specimen  subjected  to  a creep  and  recovery  test.  The 
results  Indicate  very  little  recovery;  and,  of  course, 
the  phenomenological  theory  using  strain  hardening 
predicts  no  recovery.  However,  it  should  be  noted  that 
the  nonlinear  viscoelastic  theory  would  predict  com- 
plete recovery  If  the  creep-recovery  curve  were 
extended  to  larger  values  of  time. 

figure  4 presents  the  creep  strain  response  of  a 
type  304  stainless  steel  (heat  No.  8043813)  uniaxial 
specimen  subjected  to  montonlcally  Increasing  load. 

The  results  are  In  reasonable  agreement  with  the  exper- 
imental results  and  are  very  close  to  the  results  based 
on  strain  hardening. 

figure  5 presents  the  creep  strain  response  of  a 
uniaxial  specimen  subjected  to  a reversed  cyclic 


Crtcp  Strain  ( fn/1n>  Crttp  Strain  ( Win)  _ 


Fig.  2 Creep  of  Uniaxial  Specimen 


Fig.  3 Creep-Recovery  of  Uniaxial  Specimen 


loading  history.  On  the  first  reversal,  the  nonlin- 
ear viscoelastic  theory  over  predicts  the  creep  re- 
sponse as  given  by  the  experimental  data.  This  is 
due  in  part  to  the  lack  of  significant  recovery  for 
Type  304  stainless  steel.  To  amplify  this  behavior, 
consider  the  stress  history  up  to  the  time  of  the 
second  stress  reversal, 

o(t)  - o H(t)  - 2o  H(t-t,)  , 0 < t < t2  (79) 

Substituting  equation  (79)  into  equation  (63)  gives 

tC  - 0(e)  g2  o - 2g 2 o D(*-e,).  (80) 

Thus  it  is  seen  that  the  second  term  gives  2g?  o 
times  the  compliance.  Thus  twice  the  value  ot  go  o 
is  multiplying  the  compliance  than  was  on  initial 
loading.  This  seems  reasonable  for  materials  with  an 
active  recovery,  l.e. , polymers.  However,  for  metals 
the  recovery  is  only  slight  and  the  second  term  over 
predicts  its  contribution  to  the  recovery.  Further, 
the  first  term  is  related  to  the  concept  of  a fading 
memory,  which  is  questionable  for  the  stress  reversal 
situation.  Mote  the  curve  giver,  for  the  first  stress 
reversal  with  no  recovery  or  memory  of  past  history. 
This  curve  was  obtained  by  initializing  equation  (20) 
to  zero,  which  represents  the  memory  of  the  material 
to  prior  histories,  and  setting  G11n  2 ■ 0 when  a 

stress  reversal  occurs.  This  gives  the  same  response 
as  an  initial  loading  with  the  strain  accumulated  at 
the  time  of  the  reversal  serving  as  a new  origin;  this 
is  the  current  procedure  being  followed  in  reference 
(14)  for  the  first  stress  reversal  using  toe  phenom- 
enological theory. 
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Fig.  4 Creep  Strain  Response  of  Uniaxial  Specimen  Subjected  to  Honotonically  Increasing  Load1' 
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Fig.  5 Creep  Strain  Response  of  Uniaxial  Specimen 
Subjected  to  Cyclic  Load 


CLOSURE 

The  purpose  of  this  paper  was  first  to  present 
an  efficient  method  of  integrating  Schapery's  single 
Integral  equation  for  use  In  an  Incremental  finite 
element  program  wherein  the  equations  of  motion  are 
Integrated  step-by-step.  And  cecondly  to  present 
some  preliminary  results  for  metals  using  this  theory. 

It  Is  believed  that  the  method  presented  for 
Integrating  the  single  Integral  equation  Is  efficient 
based  on  the  fact  that  the  Integrations  required  are 
for  the  current  time  step  only  and  does  not  require 
integration  over  the  entire  past  history  to  determine 
the  current  creep  rate.  From  the  experience  gained 
with  the  computer  program  NONVIS,  the  programming 
aspect  of  implementing  the  theory  Is  straight  forward 
and  can  easily  be  Incorporated  Into  existing  finite 
element  programs. 

The  results  of  this  paper  are  not  presented  to 
draw  final  conclusions  regarding  whether  Schapery's 
nonlinear  viscoelastic  theory  or  the  phenomenological 
theory  based  on  strain  hardening  Is  superior.  Obvi- 
ously, the  application  of  the  phenomenological  theory 
for  metals  has  developed  over  many  years  whereas  the 
application  of  nonlinear  viscoelastic  theory  to  metals 
Is  relatively  new.  However,  from  the  results  present- 
ed here,  it  Is  seen  that  In  all  except  the  case  of 
stress  reversal,  the  nonlinear  viscoelastic  theory 
gives  results  as  acceptable  as  those  based  on  strain 
hardening  theory  and  In  fair  agreement  with  experi- 
mental data.  For  the  stress  reversal  case,  two  prob- 
lems arise.  The  nonlinear  viscoelastic  theory  pre- 
dicts a complete  recovery  of  the  strain  accumulated 
at  the  first  and  subsequent  reversals  (as  in  linear 
viscoelasticity).  Secondly,  the  theory  has  a fading 
memory  of  past  stress  history.  However,  as  was  noted, 
these  effects  can  be  neglected  for  the  stress  rever- 
sal case  by  Initializing  to  zero  the  parameters  that 


contain  the  memory  of  past  histories  when  a stress 
reversal  occurs  and  setting  the  value  of  G,,  • 0. 

1  N-2 

Although  the  results  obtained  herein  have  been 
for  the  response  of  a stainless  steel,  most  of  the 
results  reported  to  date  have  been  for  a polymer  or 
polymer  composite.  As  noted  In  the  Introduction, 
references  (12,13)  show  the  theory  to  be  In  good 
agreement  for  a polymer  composite  for  the  following 
stress  histories:  creep,  creep- recovery , and  multi- 
step  Inputs.  Therefore,  It  Is  believed  that  the  use 
of  the  Integration  technique  for  the  single  integral 
equation  In  conjunction  with  an  Incremental  finite 
element  analysis  has  immediate  application. 

Current  work  with  the  single  Integral  equation 
of  Schapery's  Involves:  (1)  providing  a better  pre- 
dictive capability  for  the  stress  reversal  and  other 
stress  histories  for  metals,  (2)  determining  the  creep 
compliance  and  material  function  for  several  composite 
materials  and  (3)  Incorporating  the  above  Integration 
procedure  Into  a finite  element  program  for  testing  on 
two  dimensional  problems. 
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APPENOIX 

Recursive  Relationships  for  Cl^  and  C2^  ^ 

N N 

After  expanding  equation  (23)  as  the  sum  of  two 
Integrals,  It  can  be  rewritten  as 


d6cd  , cdN-l  * H-2 


■ (r)  , -<VVl>/Tr 


{/*»-*,-<  vr*'>'v  ^ 

Jo 

-(♦«.!-♦•  )/»r  dGrd  ) 

♦j  • -ardt] 


To  evaluate  equation  (81),  note  that  If  In  equation 
(23)  N Is  replaced  by  N-l,  equation  (23)  can  be 
written  as. 


■ r 

CdN-l  L 


2 -(Vl^'l/’r  d(td 


which  Is  the  value  of  Cl;v  from  the  previous  time 
step.  Cd 

Define  a*N  as, 

a*N  * *n  ’ *8-1  l 

and  approximate  the  time  rate  of  change  of  GC(J  by 


-W  («) 

on  the  time  Interval  tu.?  « t « tw.i.  Using  equation 
(82),  (83)  and  (84)  In  equation  (81)  gives, 

cl(r)  „-Vvjc,(r|  ,GcdW-i:.GcdW-2x 
CdN  l cdN-l  AtN-l 

I e " r df|  . (85) 

\-2  * 

In  the  Integral  given  by  equation  (17),  replace  N by 
N-l  to  get, 

41  • 4,  «•  ». 

*N-2 

which  Is  Jl^  from  the  previous  step.  Thus,  equation 
(85)  can  be  written  as, 

c,w.  -Vv{c,M;  0^-^w.JilJ. 

(87) 

Note  the  values  of  Clj^  , Jlirj,  G . , and  G . 

“n-1  N']  cdN-l  cdN-2 

are  known  values  from  previous  steps. 

Similarly,  equation  (24)  can  be  manipulated  and 
rewritten  as 


where  C2[r.  ^ and  J2^rJ  are  values  of  equation  (24) 
cdN-l  N*' 

and  (18)  from  the  previous  step. 

Recurvlse  Relationship  for  F^ 


After  expanding  equation  (36)  as  the  sum  of  two 
integrals,  it  can  be  written  as. 


■ss-f 


N-2  "(VVl^-l**,),lr  dGcd 


N-l  e-(V*N-l+Vr*')/Tr  dt' 


To  evaluate  equation  (89),  note  that  If  In  equation 
(36)  N is  replaced  by  N-l,  equation  (36)  can  be  re- 
written as. 
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If)  f1*’2  ,V)/l-  d6rH 

r'v/i  * -af**'  (MI 

which  Is  the  value  of  from  the  previous  step. 
Using  equation  (89),  (907  and  (84).  equation  (36)  can 
be  written  as, 

F(r)  . f-VTrfF(r)  6cdN-rGcdN-2 

CdN  i cdN-l  * 

N-l  ■(♦N.i'*,)/Tr  ) 

I « dt'J.  (91) 

\-2  * 

The  Integral  term  in  equation  (91)  can  be  written  as 
using  equation  (86)  to  obtain 

■ e’1V'r  <rf<K’ 

Evaluation  of  the  Quantities  JllrVd  J2.(r  ^ 

Tr\  — N N — 

Expand  Jl^  'In  equation  (17)  as 


'T-kiP 

^ Vl 

* e " M df  . 

Vl 

\-l 


or 

Vl  (93) 

.f" 

Vi 

Let  the  change  In  reduced  time  during  the  Interval 
^t^be  related  to  tf<  through  the  parameter,  a,  such 

*N  * Vl  • a (VVl1  or  * " Sj  • (94) 

Defining  the  quantity  Jl^r^  as 


allows  equation  (93)  to  be  rewritten  as 

^r’  ■ k,  * ■>.;« 

Vi 

■ i.  o - e'iVv)  ♦ ji:(r) . 


Making  the  following  change  of  variable  In  Jl*^r^ 

N 

T • t 't|<  + v 'Vi  , 2t''VVi 
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transforms  the  lower  and  upper  limits  of  the  Integral 
in  equation  (93)  to  -1  and  +1,  respectively,  and  gives 
dt'  = (AtN/2)dt. 
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Equation  (97)  can  be  integrated  numerically  by  stand- 
ard Gaussian  quadrature. 

Expanding  JZJL'* * ) in  equation  (18)  as 
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Upon  Integration  of  the  first  term,  equation  (98) 
becomes 
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Define, 


and  write  equation  (99)  as 
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Using  the  same  change  of  variable  as  above  allows 
equation  (100)  to  be  rewritten  as 


Evaluation  of  » 

Recalling that  the  reduced  time- Is  defined  by 

♦ •f  ^ df  ; 

Jo  0 

for  the  interval  t^  < t < tN,  one  can  then  write 
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For  small  time  Intervals,  one  can  approximate  the 
stress  as  a linear  function  of  time  with  the  Interval 
Atjj  such  that 


Chanping  the  variable  of  Integration  to  o,  equation 
(103)  becomes 
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Flg.  6 Linear  Interpolation  of  Stress 

i 

from  which 


o(t)  - A,  ♦ A,  t , t„  . < t < tM 


